Collision of Domain Walls in asymptotically Anti de Sitter spacetime 
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We study collision of two domain walls in 5-dimensional asymptotically Anti de Sitter spacetime. 
This may provide the reheating mechanism of an ekpyrotic (or cyclic) brane universe, in which two 
BPS branes collide and evolve into a hot big bang universe. We evaluate a change of scalar field 
making the domain wall and can investigate the effect of a negative cosmological term in the bulk 
to the collision process and the evolution of our universe. 

PACS numbers: 98.80.Cq 



INTRODUCTION 



The theory of inflation not only resolves the key the- 
oretical problems of a big bang theory such as the flat- 
ness and the horizon problems P, Q , but also seems to 
be conflrmed by recent observational data on CMB 
However, it is still unclear what is the inflaton in a fun- 
damental unifled theory such as string/M-theory. So far, 
there has not been found a convincing model, although 
we may have several possible links 4, _5J . 

While, a new paradigm on the early universe, the 
so-called brane world has been proposed 0, Q- Such 
speculation has been inspired by recent developments in 
string/M-theory There has been tremendous 

work on this scheme of dimensional reduction [TTl IT^ , 
where ordinary matter fields are confined on a lower- 
dimensional hypersurface, while only gravitational field 
propagates throughout all of spacetime. Among such 
dimensional reduction scenarios, Randall and Sundrum 
(RS) 12] proposed a new model where four-dimensional 
Newtonian gravity is recovered at low energies even 
without compact extra dimensions. Based on such a 
new world picture, many cosmological scenarios have 
been studie d ll^ Q, 0. See also recent reviews 
[TgL It^ ITsL I19| |. We have found some deviations from 
standard cosmology by modifications of 4-dimensional 
Einstein equations on the brane |2C1 . even for the case 
that there is a scalar field in bulk ^1]. 

In such a brane world scenario, for resolving the above- 
mentioned key theoretical problems in the big bang the- 
ory, a new idea of cosmological model has been proposed, 
which is called the ekpyrotic scenario or the cyclic uni- 
verse scenario [2^ |2^ . It could be an alternative to an 
inflationary scenario. Since it is not only motivated by 
such fundamental unified theory but also may resolve the 
above mentioned key cosmological problems, it would be 
very attractive. 
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This scenario is based on a collision of two cold branes. 
The universe starts with a cold, empty, and nearly BPS 
ground state, which contains two parallel branes at rest. 
The two branes approach each other and then collide. 
The energy is dissipated on the brane and the big bang 
universe starts. The BPS state is required in order to 
retain a supersymmetry in a low-energy 4-dimensional 
effective action. The visible and hidden branes are flat 
and are described by a Minkowski spacetime, but the 
bulk is warped along the fifth dimension. 

There has been much discussion about density pertur- 
bations to see whether this scenario is really a reliable 
scenario for the early universe [23. l25l |26L 12^ . It has 
been shown that the initial spectrum seems not to be 
produced as a scale invariant one. So this is the weak- 
ness of this scenario. However, it may turn out that we 
find other wayout to produce a scale invariant fiuctua- 
tion. Then, it is necessary to analyze the other impor- 
tant process in this scenario, e.g. the reheating process, 
in detail. For the work about the collision of branes, 
there has been some studies by |2^ with using the 
brane approximated as the delta function. In the anal- 
ysis using a delta function, the microscopic processes in 
collision such as dissipation of matter field and reheating 
of the universe can not be described. Even though there 
are some works on the reheating mechanism itself 
has not been so far investigated in detail. So the purpose 
of our work is to construct a consistent brane collision 
model and to analyze its reheating process. 

With such analysis, we may find how we can recover 
the hot big bang universe after the collision of branes. 
In our previous work (paper I) [sj, we proposed a new 
reheating mechanism in the ekpyrotic universe, such as 
a quantum creation of particles confined on the brane at 
the collision of two branes. As a thick brane, we could 
adopt a domain wall model constructed by a bulk scalar 
field and analyze the collision of two moving domain walls 
in a 5-dimensional bulk spacetime. In the paper I, we 
consider the simplest situation such that two domain 
walls collide in a given background spacetime, i.e. the 
five-dimensional (5D) Minkowski spacetime. As shown 
by several authors [32, 113, 13 HEl i the results highly de- 
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pend on the incident velocity v of the walls. We have 
shown the time evolution of colliding two domain walls 
in the paper I, which confirms the previous works. For 
a sufficiently large velocity, e.g. ■u>0.25, it is shown that 
a kink and antikink will just bounce off once, because 
there is no time to exchange the energy during the col- 
lision. For a lower velocity, we find multiple bounces 
when they collide. For example, we find the bounce oc- 
curs once for v = 0.4, while it does twice for v = 0.2. 
We also find many bounce solutions for other incident 
velocities, as shown in Appendix B of [sj (see also [s^V 
The number of the bounces sensitively depends on the 
incident velocity. A set of the values of v which give the 
same number of bounce forms a fractal structure in the 
u-space as shown Fig. 6 of [13 • Then, we evaluate a pro- 
duction rate of particles confined to the domain wall. As 
a result, the energy density of created particles is given 
as p « 20g^A^b m% where g is a coupling constant of par- 
ticles to a domain-wall scalar field, Nh is the number of 
bounces at the collision and mg, is a fundamental mass 
scale of the domain wall. It does not depend on the width 
d of the domain wall, although the typical energy scale 
of created particles is given by a; ^ 1/d. The reheat- 
ing temperature is evaluated as Tr « 0.88 g N^^^. In 
order to have the baryogenesis at the electro-weak en- 
ergy scale, the fundamental mass scale is constrained as 
m,,>l.l X 10'' GeV for g ^ 10~^. For such an energy 
scale, we may find a sufficiently hot universe after colli- 
sion of domain walls. 

In order to study whether such a reheating process 
is still efficient in more reliable cosmological models, we 
have to include the curvature effect. In particular, some 
brane universe are discussed with a negative cosmological 
constant Hence, we study here how gravitational ef- 
fects change our previous results. In order to investigate 
such an effect, we have not only to investigate the colli- 
sion of domain walls in a curved spacetime, but also to 
solve the spacetime by use of the 5D Einstein equations. 
Inspired by the RS brane model, we include a potential 
of the scalar field which provides an effective negative 
cosmological constant in a bulk spacetime. 

We first set up the initially moving two domain walls, 
each of which is obtained by boosting an exact static 
domain wall solution |^ (§11)- Although this solution is 
obtained in the four dimensions, it is easy to extend it 
to the 5D one. We then solve the 5D Einstein equations 
and the dynamical equation for a scalar field to analyze 
collision of thick walls in asymptotically AdS spacetime 
(§111). The concluding remarks follow in §IV. 

We use the unit oi c — h — 1. 



II. COLLISION OF TWO DOMAIN WALLS 

A. Basic Equations 

We study collision of two domain walls in 5D space- 
time. To construct a domain wall structure, we adopt 



a 5D real scalar field $ with an appropriate potential 
V^($), which minimum value is negative. This potential 
gives an asymptotically anti-de Sitter (AdS) spacetime 
just as the RS brane model. 

Since we discuss the collision of two parallel domain 
walls, the scalar field is assumed to depend only on a time 
coordinate t and one spatial coordinate z. The remaining 
three spatial coordinates are denoted by x. For numerical 
analysis, we use dimensionless parameters and variables, 
which are rescaled by the mass scale m$ , which is defined 
by the vacuum expectation value at a local minimum as 

■I'd = w|^', as 

$ 

t — m.s>t , z — m<;>z,^ — — jj^ . (2-1) 

In what follows, we drop the tilde in dimensionless vari- 
ables for brevity. 

It is possible to choose coordinates such that the bulk 
metric has the "2D conformal gauge" , i.e. 

ds^ = e^^(''^-'> ( ^dt^ + dz^) + e2^(*'-)da;2 . (2.2) 

This gauge choice also makes the initial setting easy when 
we construct moving domain walls by use of the Lorentz 
boost. 

In this gauge, the 5D Einstein equations and the dy- 
namical equation for a scalar field are split into three 
dynamical equations; 

A ^A" + - 35'" - k2($2 _ $'2 + ie2^F($)) , 

B =B" - 3^2 + 3B'^ + -K^e2^t/($) , 

$ =$" - 3B$ + 3B'$' - ie2^y'($) , (2.3) 

plus two constraint equations; 

BB' - A'B - AB' + 13' = , 

3 

2B'' + B" - A'B' -AB-B^ 

= --kI{4>^ + <i>'^ + e'^^V{<P)), (2.4) 
3 

where a dot (') and a prime (') denote d/dt and d/dz, 
respectively. 

These are our basic equations. Before solve them nu- 
merically, we have to set up our initial data, which sat- 
isfies the constraint equations (|2.4|l . 

B. Domain Wall Solution 

For an initial configuration of a domain wall, we use an 
exact static solution given by [36j . They assume a scalar 
field <& with a potential 
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where 



with 



W 



(2.6) 



is a superpotential, and k\ and d are the five dimen- 
sional gravitational constant and the thickness of a do- 
main wall, respectively. The potential minima are lo- 
cated at $ = ±1 in the range of |<1>|<5 (see Fig. 01. The 
potential shape is similar to a double-well potential, but 
it is asymmetric. The minimum value at $ = 1 vanishes, 
while that at $ = —1 is negative. With this potential, we 
can obtain analytically an exact solution for two colliding 
domain walls as follows. 




FIG. 1: The scalar field potential is plotted where 

d = \/2, K5 — 0.3. This potential behaves like a double-well 
potential for |<J>| < 5, except that one of its two minima has a 
negative value. On the other hand, for |$| > 5, this potential 
behaves differently from a double-well potential, that is V{^) 
rapidly falls into — oo. 



First, we show a static domain wall solution with this 
potential. A kink solution of a scalar field {K) is de- 
scribed as 

$K(y)=tanh(^) , (2.7) 

and the metric of 5D spacetime is 

ds^ = e^'^''^y\-dt^ + dx^) + dy^ , (2.8) 

with 

tanh^(y/d) 



AK{y) = --4{\n [cosh(^) 



-I- 

di 

(2.9) 

Since this exact solution is not given in our gauge, we 
have used a new coordinate y, which will be associated 
with z later. 

This metric approaches that of the AdS spacetime in 
one asymptotic region (y <C —1), i-e. 



with 



e-^/^lyl as y 



9d 



(2.10) 



(2.11) 



While it becomes a flat Minkowski space in another 
asymptotic region (y ^ 1), i.e. 

„2Ak , „2A„ 



l4 (log2-i 



(2.13) 



We depict the behaviour of metric function exp[Aif(?/)] 
in Fig. El 
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FIG. 2: The metric component of the exact solution for a 
static domain wall is plotted where we set d — \/2 and 
K5 — 0.3. This spacetime approaches Minkowski space be- 
cause Ak Aoo (a constant), while it becomes asymp- 
totically AdS because Ak ^k\y\ as y ^ — oo, where 
k = 8/ti/9d. 



By reflecting the spatial coordinate y, we also flnd an 
antikink solution {K) as $/f(y) = ^K{—y) ~ —^K{y)- 
The corresponding metric of this antikink solution is also 
obtained by reflection of y-coordinate, i.e. Ai^{y) — 
AA-(-y) = -AK{y). 

In order to describe this solution under our gauge con- 
dition l|2.2() . that is, in the {t, z) frame, we should trans- 
form the present (i,y) coordinates fEa. (|2.8|) ') into the 
(t, z) ones fEa. (|2.2fl ') by defining the coordinate z as 



dy. 



(2.14) 



This integration will be performed numerically to find 
initial data of collision of two domain walls. 



C. Moving Domain Wall 

When a domain wall moves with constant speed v in 
the fifth direction z, we can obtain the corresponding 
solution by boosting a static kink solution [K) as 



$^,(z, t) — tanh 



1 



y* {j{z - vt)) 



(2.15) 



as y 



oo . 



(2.12) 



where y* and z* are comoving coordinates of a domain 
wall, and y* (z*) is obtained by the inverse transformation 
of Eq. H2.14|l . The Lorentz transformation gives z* = 
7(z — vt) where 7=1/ y/1 — v'^ is the Lorentz factor. We 
have assumed that the center of a domain wall is initially 
located at z = 0. 
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The corresponding metric is easily obtained by Lorentz 
boost. Because of the Lorentz invariance in our 2D con- 
formal gauge, i.e., —dt*"^ + dz*"^ — —dt^ + dz'^, we find 

dsljy = exp[2AK (7(2 - vt))]{-dt'^ + dz^) , (2.16) 

where Ax{z*) = A^iy* (z*)). The function Ax{y) is 
given by Eq. H2.9|l . The center of a domain wall (z* — 0) 
moves as z = vt in our {t, z)-coordinate frame. Then 
we regard that the metric describes a spacetime with a 
domain wall moving with constant speed v in the z di- 
rection as well as a scalar field $ does so. 

In order to discuss collision of two moving domain 
walls, we first have to set up its initial data. Using Eqs. 
(|2.15l) and (|2.16|) . we can construct an initial data for 
two moving domain walls as follows. Provide a kink so- 
lution a,t z = —zq and an antikink solution at z = zq, 
which are separated by a large distance and approaching 
each other with the same speed v. We then obtain the 
following initial data; 

0) - + zo, 0) - -zo,0)-l, 
A{z,0) = A„(z + zo,0)- A_^(z-zo,0)- Aoo, (2.17) 

where Aoo is the constant value given by Eq. H2.13|l . The 
initial values of $ and A are also given by 

$(z,0) = ^^{z + zo,0) - zo,0), 

A{z, 0) = i„(z + zo, 0) - i_„(z - Zo, 0) . (2.18) 

Obviously, we set A = B and A = B at initial. 

The spatial separation between two walls is given by 
2zo, and as long as the separation distance is much larger 
than the thickness of the wall {zq ^ d), the initial con- 
ditions ifTTTjl and (|TTH|) give a good approximation for 
two moving domain walls. Using these initial values, we 
solve the dynamical equation H2.3|l numerically. The re- 
sults will be shown in the next section. 



III. TIME EVOLUTION OF COLLIDING 
DOMAIN WALLS 

We use a numerical approach to solve the dynamical 
equations (|2.3|) for the colliding domain walls. We adopt 
a numerical method similar to one used in [3]|. The dif- 
ference is found in a boundary conditions. We impose the 
Dirichlet boundary condition for the scalar field, which 
is the same as the paper I, while the Neumann bound- 
ary condition is used for the metric as A'{z) = —k^e'^^^\ 
which is derived from the asymptotic form of the metric, 
i.e., e^(^) ^ l/(fc(7l^| + 1)) as \z\ 00. 

We have three free parameters in our model, i.e. a 
wall thickness d, an initial wall velocity u, and a warp 
factor k (or the gravitational constant K5). Two of them 
(d, fc(or K5)) are fundamental because they appear in the 
theory. In the paper I, we studied the collision of two 



domain walls in the fixed Minkowski background |32l |. 
where we had two free parameters d and v. So, including 
the gravitational back reaction, we investigate how K5 (or 
k) changes the previous results. In what follows, fixing 
the value of d, i.e. d = -^72, we show our results. 



A. Time Evolution of Scalar Field $ 

First let us set K5 = (or k = 0), that is the case of 
Minkowski background spacetime. Although the scalar 
field potential is slightly different from that in the paper 
I, the result is exactly the same. This simulation also 
gives a check of our numerical code. 

Next we perform our simulation for K5 7^ 0. For a small 
value of K5, i.e., K5<0.05 (equivalently fc<1.57 x 10"'^ 

or m$<(0.05)^/^m5 « O.lSBms), the collision process is 
very similar to the case of the Minkowski background. 
Setting the initial velocity v — 0.4, we show our numeri- 
cal results for K5 ~ 0.05 in Figs. 01 and 01 The evolution 
of $ is depicted in Fig. 01 while that of the energy density 
is shown in Fig. ^ 




FIG. 3: Collision of two domain walls where the initial veloc- 
ity V = 0.4. The time evolution of the scalar field $ is shown 
from t = to 70. The collision occurs once around t = 31. 
We set d = %/2, K5 = 0.05. This process is similar to the 
Minkowski case, K5 = 0. 



The energy density of the scalar field is given by 

p^ = e-^^(^^^ + +V{'P). (3.1) 

We find some peaks in the energy density, by which we 
define the positions of moving walls (z = ±zw(t)). If 
a domain wall is symmetric, its position is defined by 
$(2:) — 0. However, in more general situation, just as in 
the present case that the scalar field is oscillating around 
some moving point, it may be natural to define the posi- 
tion of a domain wall by the maximum point of its energy 
density. 
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FIG. 4: The time evolution of scalar field energy density in 
the case of Fig. |21 The maximum point of defines the 
position of a wall {z = z^f/{t)). 
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FIG. 5; Time evolution of the position of one brane z = z-w 
which starts to move at z = — 13 with a constant speed v = 
0.4. We set d = V2, ks = 0.05. 



Fig. |31 denotes the position of brane z = z^f^it) with 
respect to t. The brane moves with constant speed v = 
0.4 toward the colhsion point z = 0, and coUide, then 
recede to the boundary. We also find smaU oscillation 
around a uniform motion after collision. 

Since we assume that we are living on one domain wall, 
we are interested in a particle '0 confined on the domain 
wall. If a particle -0 is coupled with a 5-D scalar field 
$, which is responsible for the domain wall, we expect 
quantum production of V'-particles at collision of domain 
walls. This is because the value of the scalar field $ on 
the domain wall changes with time. This fact may play 
an important role in a reheating mechanism 31]. Once 



we find the solution of colliding domain walls, we know 
the time evolution of a scalar field on the domain wall, 
and we can evaluate production rate. 

At the position of a domain wall, the induced metric 
is given as 

ds^ =-dT'^ +a'^{T)dx'^ , (3.2) 
where proper time is determined as 

r = j e^^dt, (3.3) 

and the scale factor is 

a(T)=e^*(T). (3.4) 

Aw and i?-vv are evaluated on the brane, i.e., Aw ~ 
A{t, zyj{t)) and i?w = B{t, zw(i)). The Hubble parame- 
ter of the brane universe is defined by 

H{r)^l^ ^^^e-^-BMr). (3.5) 

a ar z=zw dr 
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FIG. 6: Time evolution of a scalar field on the moving wall 
with V = 0.4, d = V2. We set K5 = 0,0.05,0.1. The value 
of the scalar field on the wall is defined by $w ~ ^{t, 2w(i)), 
where zw{t) is the position of the wall. The time of period 
of oscillation gets slightly longer and the amplitude is a little 
bit larger as K5 increases. 

For V = 0.4, we depict the time evolution of $w on one 
moving wall for different values of K5 in Figs. IHIandEI 
The feature of collision is similar, but the behaviour of a 
scalar field on the moving wall after collision is different 
for each K5. 

We summarize our numerical results for each value of 
K5 in order. 
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(i) K5=0.01 

For a small value of K5, e.g., K5 = 0.01 (equivalently 
k ^ 6.29 X 10~^ or to$ ^ 0.0464to5), the result is al- 
most the same as the case of the Minkowski background, 
in which case we find one bounce point, which corre- 
sponds to a crossing point in Fig. 01 and then the os- 
cillations around — follow (see the dotted line in 
Fig. This oscillation is explained by using a pertur- 
bation analysis in Minkowski spacetime 31]. We have 
found one stable oscillation mode around the kink solu- 
tion. This oscillation appears by excitation of the system 
at collision. 

(h) K5 = 0.05 

As increasing the value of K5 slightly larger, for ex- 
ample K5 = 0.05, (equivalently for k ^ 1.57 x 10"'^ 
or TO$ ~ O.I36TO5), the time evolution of $w slightly 
changes. In Fig. El just as the case of K5 = 0, we find 
one bounce and successive oscillations. However, the pe- 
riod of oscillation is slightly longer and the amplitude 
gets a little bit larger as K5 increases. 

(iii) K5 = 0.15 

For K5 = 0.15, (equivalently k ^ 1.41 x 10~^ or m$ ~ 
0.282m5), the behaviour of this oscillation changes dras- 
tically. After several oscillations, the scalar field leaves 
= as shown in Fig. [7| The numerical simulation 
eventually breaks down because all variables diverge. 
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FIG. 7: Time evolution of a scalar field on the moving wall 
for V = 0.4, d = ^/2, K5 — 0.15. The value of the scalar field 
is given by <l>w = ^{t, zw{t)), where zw{t) is the position of 
the wall. "l>w goes out of oscillation phase at r ~ 55 and 
then numerical simulation stops. 



(iv) K5 > 0.15 

For the larger values of K5 than 0.15, the time to appear- 


















\ 






i 











(a) K5 = 0.2 



(b) tZ3 = 0.25 



FIG. 8; Time evolution of the scalar field $w for v — 0.4, d — 
\/2, setting (a) K5 = 0.2 and (b) K5 = 0.25. The numerical 
simulation breaks down when (a) t ~ 42 and (b) r ~ 35, 
respectively. 



ance of singularity becomes shorter, that is, contrary to 
the case of K5<0.15, the scalar field after collision does 
not oscillate but leave $w = soon. The time evolution 
of $w is shown in Fig. H for K5 = 0.2 {k = 2.51 x 10-^) 
and K5 = 0.25 {k = 3.93 x 10"^). 

The metric component ^ at z = also diverges as 
shown in Fig It is not a coordinate singularity, but 
a curvature singularity. In order to show it, we calcu- 
late the so-called Kretschmann invariant scalar, which is 
the simplest scalar invariant quadratic in the Riemann 
tensor, and is defined as 



T^abcd r> 
n nahcd 



-4A 



3{B + - AB- A'B' f 
'? 

3{B'" - B^f + 'i{B" + B'" -AB- A'B'f] . (3.6) 



3(yl'B + AB' -B' - BB'f + {A - A" f 

'2 D2^2 



In Fig. ini^b), we depict the time evolution of the 
Kretschmann scalar at the origin z = 0, which diverges 
at t ~ 69. It is caused by the divergence of a quantity A. 
We conclude a singularity forms at the origin z = 0. 

This divergence is not a numerical error because a con- 
straints equations (|2.4() are always satisfied within 10~^ 
- 10~^ % accuracy except at time when the singularity 
appears. 

In Appendix 1X1 we study the reason why a spacetime 
is unstable and eventually evolves into a singularity in 
detail using a perturbation analysis of the Einstein equa- 
tions and dynamical equation of <&. We find that a per- 
turbed oscillation mode around an unperturbed kink so- 
lution becomes overstable for K5 = 0.1. From our anal- 
ysis, we conclude that gravitational back reaction makes 
a kink solution unstable contrary to the Minkowski case. 

Next, we show the result for the case of the initial 
velocity w = 0.2. In the Minkowski case, this incident 
velocity shows two-bounce at collision process (see Fig. 
El (a)). After two walls collide, they bounce, recede to 
some finite distance, turn back and then collide again. 
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FIG. 9: (a) Time evolution of the metric component A at 
z = for 11 = 0.4, d = ^2, K5 = 0.15. It diverges at f ~ 70. 
(b) Kretschmann scalar invariant {R°'^'^'^ Rabcd) at the origin 
z = Q for the case of f = 0.4, d — \/2, K5 — 0.15. We find that 
the Kretschmann scalar diverges at t ~ 69. This means that 
it is not a coordinate singularity, but a curvature singularity. 
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(a) K5 = (Minkowski) 



(b) K5 = 0.05 



FIG. 10: Time evolution of a scalar field $w for v — 0.2, d = 
y/2. For (a) K5 = and (b) K5 = 0.05, we find two peaks which 
correspond to twice bounces at collision. Moreover, it is seen 
that an effective negative cosmological constant prolongs the 
time interval between two bounces. 



For small values of K5, e.g., K5 < 0.05, the collision pro- 
cess is very similar to the case of K5 = (see Fig. ITUT b')'). 
As K5 increases, the time interval between first and sec- 
ond bounces becomes longer as shown in Fig. This 
can be understood from the fact that the above men- 
tioned oscillation after collision will radiate the energy. 
So a kink-antikink pair is loosely bounded and it takes 
longer time to collide again. 

For the case of K5>0.1, this feature of collision is dras- 
tically changed. Two-bounce collision never occurs, that 
is, two walls collide only once as shown in Fig. E| This 
is because a lot of energy of a kink-antikink pair is radi- 
ated away via the unstable oscillation after collision and 
it has not enough energy to form a trapped state. After 
the first bounce, the domain walls never collide again but 
recede each other. 

For larger value of ^5, we find only one-bounce col- 
lision. Namely, a "negative cosmological constant" out- 
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FIG. 11: Time evolution of the position of one brane z = 
zw which starts to move at 2: = —13 with a constant speed 
V — 0.2, setting d — \/2 K5 — 0.05. The dashed line denotes 
Ks = 0.05 while the dotted line does K5 — 0. The time interval 
between two bounces is prolonged by an effective negative 
cosmological constant. 



side a kink-antikink pair keeps away two walls toward the 
boundary, so it plays as an effective attractive force. 



B. Time Evolution of Metric 

We evaluate the time evolutions of the metric Aw, Bw 
on the brane and plot them in Figs. El Both of two 
quantities decrease with time except that _Bw increases 
slightly through the bounce. 

From those two quantities, using Eqs. H3.4|) and (|3.5|l . 



we evaluate a scale factor of our universe a(T) 



and the Hubble expansion parameter H = ^/a, where 
T is the proper time of domain wall defined by Eq. , 
and show them in Figs. El setting v = 0.4. From this 
figure, we find that our universe expands slightly before 
bounce then eventually contracts. For each K5, the scale 
factor a and the Hubble parameter H are plotted in Figs. 
1151 and [TCI From these figures, we see that our universe 
contracts faster as K5 gets larger, i.e. a negative cosmo- 
logical constant increases. 

Next we investigate the scale factor a and the Hubble 
parameter H for the case oi v = 0.2. Setting K5 = 0.05, 
that is the case of two-bounce collision, a and H are 
plotted in Fig. El We find two contracting phases, which 
correspond to each bounce at collision. 

For K5>0.1, the bounce occurs once only because of a 
large negative cosmological constant, a and H are plot- 
ted in Fig. El for K5 — 0.1. In this figure, we find the 
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FIG. 12; Time evolution of the scalar field $w for v = 0.2, 
d — \/2, = 0.1. Two-bounce process does not occur. 
Compare it with the cases in Fig. 1101 



FIG. 15: Time evolution of the scale factor a = e^™ for 
V = 0.4, d = V2. We set K5 = 0.1,0.15 and 0.2. As K5 gets 
larger, the speed of contraction becomes larger. 




FIG. 13: Time evolution of the metric A and B on the moving 
wall for V = 0.4, d = \/2, 1^5 ~ 0.15. The value of the metric 
is given by Aw ~ A{t, zwit)), where zw{t) is the position of 
the wall. Both of two quantities decreases with time. 




(a) a 



{h)H 



FIG. 14: Time evolution of the scale factor a = e ™ and the 
Hubble parameter H = d{T)/a for v — 0.4, d = \/2, ^5 — 
0.15 with respect to the proper time r. 
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FIG. 16: Time evolution of the Hubble parameter H = a{T)/a 
for V = 0.4, d = a/2. We set K.5 = 0.1, 0.15 and 0.2. As K5 gets 
larger, the typical time scales of expansion and contraction 
become larger. 



universe contracts slower than the case v — 0.4. Our 
universe contracts faster as v gets larger. 

Finally, we find in Fig. ^]that there are two discon- 
tinuous stages in the evolution of the Hubble parameter 
at t ~ 30 and ~ 32. As we will see, we can conclude that 
these discontinuities appear just because of ambiguity of 



9 




(a) 



{h)H 



FIG. 17: Time evolution ol the scalar factor a and the Hubble 
parameter H on the moving wall for v = 0.2, d — \/2, K5 — 
0.05. We find two contracting phases. 




(a) a 



{h)H 



FIG. 18: Time evolution of the scalar factor a and the Hubble 
parameter H on the moving wall for v = 0.2, d = y/2, K5 — 
0.1. The dashed and dotted lines denote the cases of v — 0.2 
and of V — 0.4, respectively. We find the universe contracts 
slower in the case of t> = 0.2 than in the case v = 0.4. 



the definition of a wall-position, z = z-yv. The Hubble 
parameter mainly depends on the derivative of the met- 
ric with respect to t, B. This quantity in fact has two 
discontinuous stages as seen in Fig. 1191 In this paper, we 
define the position of wall by one where the energy den- 
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FIG. 19; Time evolution of B and iw for v = 0.4, d — 
\/2, K5 — 0.15. We find two discontinuous stages at t ~ 30 
and t ~ 32. 



sity of scalar field gets maximum. However, this position 
does not move continuously through a bounce. Actually, 
we show iw has also two discontinuous stages as shown 
in Fig. ^1 where the speed of a wall apparently exceeds 
the speed of light. Namely near the bounce, the defini- 
tion of wall-position is not well-defined. This is because 
there is no wall configuration during the collision. Hence 
these discontinuities of the Hubble parameter seem to be 
apparent. We should look at the global time evolution. 



IV. SUMMARY AND DISCUSSION 

We have studied collision of two domain walls in 5D 
asymptotically Anti de Sitter spacetime. This may pro- 
vide the reheating mechanism of an ekpyrotic (or cyclic) 
brane universe, in which two BPS branes collide and 
evolve into a hot big bang universe. We evaluate the 
values of both a scalar field corresponding to a domain 
wall and metric on the moving wall for different value 
of the warp factor k which is related to a gravitational 
effect K5. We analyze two typical incident velocities, i.e. 
V = 0.4, and v — 0.2, which correspond to one-bounce 
and two-bounce solutions in the Minkowski spacetime, 
respectively. 

For the case oi v = 0.4, the global feature of collision 
does not change so much for different values of K5, but 
the behaviour of oscillation after the collision is different 
for each K5. For small value of K5<0.01, the oscillation is 
the same as Minkowski case, but for K5>0.05, it becomes 
an overstable oscillation. So its period and amplitude 
get larger as K5 increases. In the cause of this unstable 
oscillation, the singularity appears after collision. This 
singularity is very similar to that found in Khan and Pen- 
rose (37j , in which they discuss collision of plane waves 
and formation of a singularity. Hence the appearance 
of singularity in the present model could be understand- 
able because we take into account a gravitational effect 
in collision of two domain walls. 

In the time evolution of our universe, we find that the 
universe first expands a little just before collision and 
then contracts just after collision. This result is consis- 
tent with [25. We cannot explain our hot big bang uni- 
verse as it is. It is also found that the speed of expansion 
and contraction gets faster as K5 increases. 

For the second case, i.e., v — 0.2, we show the bounce 
does not occur twice for larger value of K5 (k5>0.1) cor- 
responding to the unstable oscillation. 

We shall discuss about the value of a warp factor k. 
We consider a curvature length I — \/k written in the 
following form 



/ = 1.97 X 10 



10- 



-)( 



TeV \ 



m 



(4.1) 



where m$ is a mass scale of a domain wall. Here we set 
a value of a warp factor k in the region 0.01<fc<0.25. 
On the other hand, we also know from the experimental 
data of testing a gravitational inverse-square law that the 
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curvature length must be smaller than 0.1 mm '38]. From 
this constraint equation, we obtain 



fc>1.97xl0-i5f— ) 



(4.2) 



So the values of k used in our simulation satisfy this 
constraint. 
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APPENDIX A: PERTURBATIONS OF A 
DOMAIN WALL 

In Fig. El we find oscillations after collision and those 
amplitudes and periods increase as K5 gets larger. To 
understand this feature, we analyze perturbations around 
a static domain wall solution in this appendix. 

We use the coordinate y, by which a static domain wall 
solution is given by analytically {^K{y), AK{y)), which 
are given by Eqs. (|2.7|) and 1)2. 9|l . We perturb the basic 
equations H2.3|l and (|2.4() by setting A — (y) -I- a{t, y), 
B = Axiy) + 6(t, y), and $ = $i^(y) -f cj){t, y). We find 
two sets of perturbation equations: 
(1) dynamical equations 



-2Au 



-2Au 



-2Af 



b = 



CI a 

dy"^ 

dy"^ 
1 

^2 



dAx da (IAk db 



dy dy 

dy dy 
dAx db 



6- 



dy dy 



ldV_ 
3d$ 



7 



dy dy 3 



2V\^a- 



K 

dV 



,(A1) 



K 



.dA 



K 



d^K db 



dy dy dy dy 



2'^ 



d^V 



K 



K 



(A2) 



(A3) 



(2) constraint equations 

db dAx . 2 



dy 



dy 

, dAK db 



d^ 



K ; 



dy 
dAK da 



dy dy dy dy 
V\^a-\ 6 + 



K 



dy dy 



(A4) 



(A5) 



In order to find the eigenvalue and eigen functions, we 
set a = a(?/)e'"*, b = 6(y)e*'^*, and = 0(?/)e*"*. Then 
the constraint equation (|A4I) is reduced to be 



db 
dy 



dAK . 
dy 



2 ^d^K ~ 



dy 



(A6) 



Inserting Eq. l|A6p into another constraint (|A5|) and us- 
ing the equations for a background solution, we find that 
the constraint equation (|A5I) turns out to be trivial. So 
we have only one constraint equation l|A6|l . 

Eliminating db/dy in HAljl and (|A3|I by use of Eq. I|A6)I . 
we obtain two coupled perturbation equations in terms 
of a, (f) as 



fd 
dy^ 



dy'^ 



[dA 



K 



\ dy 

dAK d^K 
dy dy 



-2A1, 



K 



1 dV 

3 d¥ 



K 



2Ki 



,d<i> 



K 



dy dy 



dAK da 
dy dy 

(A7) 



dy_ 



K 



dAK 
dy dy 

Id^V 



dA 



K 



dy dy 



2 



-2Ai 



UJ 



(A8) 



Eq. (|A2p is guaranteed by the other two dynamical equa- 
tions and constraint equations. Eqs. HA7|I and (|A8|I have 
the asymptotic forms as y ^ cxd as 



where 



Ai 
A2 



Id'^V 



2 d$2 



$=1 



(A9) 
(AlO) 



(All) 
(A12) 



Here we choose both negative signs in Eqs. (|A9() and 
(|A10p because negative signs correspond to out-going 
wave modes. 

We solve numerically these equations l|A7l IA8p con- 
necting the above asymptotic solutions and find the com- 
plex eigen frequency uj. For K5 — 0.1, we obtain a stable 
mode as oj^ — 1-23 -f 1.07 x and unstable mode as 

Wu = 0.644—2.90 X 10^^^, Compared this unstable mode 
with the value obtained from the oscillations after colli- 
sion found in Fig. (Notice that we use a proper time 
T in Fig. El Then it should be evaluated in the physi- 
cal time t). We obtain 0.772 « 1.2 x WyOJ^] for the real 
part of the frequency and —0.029 ~ 1.0 x ^[wu] for the 
imaginary part. So we may conclude that the overstable 
oscillations after the collision of two domain walls found 
in Fig. El are explained by the unstable mode around a 
static domain wall solution. 
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